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One of the achievements of operations research is to help managers make the right decisions. In the 
meantime, one of the operations research techniques that helps managers to achieve their goal in the 
optimal form (maximum or minimum) of a linear function under constraints of linear equations or 
inequalities is the technique of Linear Programming (LP) in definite conditions [1], [2]. In a typical 
LP model, all input parameters of the model are known values (fixed and non-probable). In real cases, 
this assumption rarely applies. The LP model is usually formulated to select some future activities. As 
a result, the parameters used will be based on the prediction of future conditions, which will inevitably 
include degrees of uncertainty [3], [4]. Due to the complexities and scarcity of information in the real 
world, in recent decades, the design of methods for recognizing, modeling and managing uncertainty 
in the behavior of systems has become one of the fascinating scientific topics and various theories 
and methods such as fuzzy set theory and probability theory for The study of uncertain systems has 
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simplex algorithm 


been developed [5]. In fuzzy set theory, experts ate used and in probability theory, sampling is used. 
However, if the number of experts and the level of experience are so low that it is not possible to extract 
membership functions or the number of samples is low, these theories can no longer be used. To address 
such issues, Professor Joe Lang Deng first proposed a different model in 1982 called the Grey System 
Theory (GST) [3]. Grey Linear Programming (GLP) is a model of analysis of Grey Systems (GSs) for 
decision making in conditions of uncertainty, which is an extension of the conventional LP method [6]. 
The following can be considered of the methods proposed by researchers to solve the problem of LP with 
grey parameters. The first category is methods that use the concept of Grey Number (GN) whitening to 
solve LP problems with grey parameters. In this method, the grey information of the GLP problem is 
written as white or definite (real) numbers using the GN whitening methods and the GLP problem is 
written as the ordinary LP problem. In these methods, the answer to the initial problem is obtained by 
solving several ordinary LP problems. Among these methods, we can mention [6]-[12]. These methods can 
be used for any kind of problem and real number operators. However, some input information may be 
lost in the parameter whitening process and the resulting answer may not be the original answer. Also, in 
these methods, the final answer is not necessarily grey and the uncertainty in the input data is not well 
reflected in the output answers. The second category is the method of finding the answer to the GLP 
problem based on the use of overlapping numbers and inverse concepts of the grey matrix. This method 
is generally a good idea for solving GLP problems [13], [14]. Excessive calculations and failure to achieve 
the stop condition after finding the answer are the two main drawbacks of this method. The third category 
is amethod based on the concepts of grey prediction [15]. In this method, we first obtain the desired values 
using the grey prediction method and then solve the problem of ordinary LP obtained. The fourth category 
is a method that uses the display of GNs in the form of intervals and the ranking of intervals to solve GLP 
problems. In this method, by converting the objective function of the GLP problem into several objective 
functions and defining the constraints related to each objective function, we separately turn the problem 
into several ordinary LP problems and solve each one [10], [16]. These methods are able to solve GLP in 
different conditions. Different sequence relationships are used to compare GNs to determine the input or 
output variable. Some researchers use the concepts of LP solution of interval numbers [6], [17]. The 
problem with these methods is that they ignore the grey space of the problem. In these methods, they 
solve the problem by using the upper and lower boundaries of the answer space, which in some cases may 
not have one of the problems. The fifth category is a method in which an attempt is made to solve the 
GLP problem directly and without the need to bleach the parameters of the problem using the simplex 
method, thus reflecting the uncertainty in the input data better in the final answer [18]. In this method, for 
comparing GNs, simultaneous comparison of center and degree of grey is used, which has a better 
performance in distinguishing between GNs. The disadvantage of this method is that it is used only to 
solve GLP problems with the grey objective function. In this paper, an algorithm such as the initial simplex 
algorithm is proposed to solve interval GLP problems, for which we do not need to convert the GLP 
problem to the classical LP problem. 


The article is presented in 5 sections. After the introduction, in Section 2, the theory of the GS and the 
concepts related to interval GNs and their comparison are stated. In Section 3, the problem of GLP 
problems is introduced. In Section 4, the simplex method for GLP problems (with objective function 
coefficients and grey sources) and the proposed algorithm are presented. In Section 5, we present and solve 
a numerical sample to express the effectiveness and understanding of the proposed method. Finally, 


Section 6 contains the conclusion. 


2 | GSs Theory 


Various theories have been proposed to deal with the uncertainties in various subjects of real-life everyday 
life, including fuzzy system theory, probability theory, and so on. GSs theory is one of the methods 
proposed to study the uncertainty of problems that are highly inaccurate due to limited data and limited 
information. The GS is defined as a system that contains uncertain information [19]. GS theory has 
attracted many researchers [20]-[29]. 


In this section, the necessary definitions for studying the GS are considered. 


Definition 1. The interval GN is defined as follows [30]. 
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@xe[x,x]={tx<t<x}, x< x. (1) y) 


Definition 2. [31] Let @x, €/x,,x,/ and @x, €/x,,x,/ be two interval GNs. 


O27 


I. For any interval GN, the center @X is defined as; 


X+X 
@x==_. (2) 
2 
II. For any interval GN, the length 7 (@x) is defined as; 

¢(@x) =|x-x|. 3) 3 
S 
IH. For any interval GN, the degtee of greyness g* (@x) is defined as; a 
A 
£(@x & 
g (®x)= ( ere (4) : 
@x 2 
ig 
IV. The main arithmetic operations can be defined on GNs. E 
fs) 
@x, + ®x =| x, +x,/X, +, | 5) F 
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2.1 | Comparison of GNs 
Comparing GNs with each other is very important for making the right decision in a grey environment. 


Dervishi et al. [21] have compared GNs in more detail. In this article has been used the concept of 
kernel and degree of grey to compare interval GNs. 


Definition 3. [32] For two interval GNs @x, and ®x,: 


@x, < @x, > @x, <c @x,. (9) 
e (x, ) =f (x, ) => ®x, =, @x, 
@BK, = OX, > 4 g°(@x,)<g"(@x,) > Ox, >, Ox, . (10) 
oa (x, ) >s (x, ) => @x, <, @x, 
3 | GSs Theory 


The problem with GLP is generally as follows. 


Maximize ®z=, )/®c, ®x, 
Subject to ™ 


n 


@a, @x, <, @b,, i=1,2,...,m, (11) 
83 " @x,2,@e, j=1,2,...m, 
@c,,@a,, @x,, @b, eR(8), i=1,2,...,.m,j=1,2,...,n. 


this section, we define LP problems involving GNs as follows: 


Maximize ®z=, ) @c, @x, 


Subject to 
n (12) 
z pa a, ®x, <®b,, i=1,2,...m, 
2 - @x,20, j=1,2,...0. 
FS 
Ey 
3 Its matrix form will be as follows: 
é Maximize ®z =, ®C ®X 
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2 13 
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Definition 4. Any vector ®x of interval GNs that satisfy the constraints of the GLP in Eg. (72) is called 
a feasible solution. Assume that Q the set of all feasible solutions to the GLP in Eg. (12). Then ®x, ¢ Q 


is said to be an optimal solvable solution to the GLP, if for all®@xeQ; @C @x<, @C @x,. 


3.1 | Basic Feasible Solution 

Consider the GS in Eg. (13), where A is a mxn matrix, ®b and @x are an m vector and n vector, 
respectively, satisfying rank( A, @b) = rank( A) =m. After possibly rearranging the columns of A, let 
A=[B,N], whete Bis a mxm invertible matrix and NV is a mx(n-m) matrix. The solution 


OX =| @x},,Ox4, | to the equations A @ X=, @b, where @X, =x, Ox, -... Ox, | and 


@X,, =¢ @0, is called a basic solution of the system. If @x, 2. @0, then ®X is called a basic feasible 
solution of the system and the corresponding grey objective value is @Z=, @C, ®X,, where 


@C, = Oy p- OCy | - Now, corresponding to every index I<j<n, define: y,= Bia, and 


@z, =, @cy Y;- Observe that for any basic index j=B,1<j<m, we have: 


1 
@z,- Be, => @c, B a,—@e, => @c, e,— 8c, => @c,-@e, => Wo. 


Where e, is the ith unit vector. Note that B is called the basic matrix and N is called the non-basic 


matrix. The components of ®X;, are called basic variables and the components of ®X,, are called non- 


basic variables. 


Theorem 1. [24] If there is a basic feasible solution with grey objective value @Z such that 
Oy =ee4, —8C, <, @0 of OZ <, @C; for some non-basic variable @x,, and ¥,=B ‘a, £O 
1< j <n, then it is possible to obtain a new basic feasible solution with a new grey objective value @Z' 


that satisfies @Z <,, @Z’. 
Theorem 2. [24] If there is a basic feasible solution satisfying 
QV, = @C, B! a, -@C, =, @Z, -@C, <, @0, for some non-basic  vartiable®@®x,, and 


Va kc 0, 1=12,..,mM, then the problem in Eg. (12) has an unbounded optimal solution. 


ik ~G 


Theorem 3. The basic feasible solution @X, =, B’@®b, @X,, =, @0 is an optimal solution to the Eg. 
(12), whenever @c, B’ a, 2,00, forall j= 7, 27. 
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simplex algorithm 


Proof. Suppose @X7 =[@x},@x4, Ja basic feasible solution is for the GLP Eg. (12). So that 


@X,=, B’@b and @X,, =, ®e, therefore, the optimal value of the objective function of Eg. (12) for 


the solution @X? will be as follows: 
@z, =, @C, @x! =, @C, ®X, =, @C, B’ @b. 


On the other hand, for each basic feasible solution @x , we have the following relation: 


A @X<, @b=, A@X=, B@X, +N @X,. 


Therefor; 


@z =, @C @X =, @C, @X,+@C, @X, =, @C, B' @b- 
> (@C, Bt a,-@C,)®X,. 


J#B; 
And, 
z=, ®z, — }}(@z,- @C, JX, 


j#B; 


According to this relation and Theorem 1, the proof is complete. 
4 | Grey Linear Programing Problem Based on the Primal Simplex 


Consider the GLP problem as in Eg. (12). Suppose the matrix B is base, in this case we rewrite the GLP 
problem as: 


Maximize ® z => @C, @X,+@C, @X, 
Subject to 
B@X,+N®X, =, ®b, 
@X, Za @X, Ear 


Hence we will have: 
@X, +B" N@X, =, B’@b>®@X, =, B'@b-B" N®@X,. 
Therefor: 


@Z=, @C,| B'@b-B" Nox, recy @X, 
=, ®C,B1®b-@C, B1 N@®X, 4+ @C, @X, 


@Z sfec, BI NaC, @X, =, @C, B1®b. 


Table 1. Simplex table related to GLP problem 


Basic Variables @X, @X,, Sources 


o : @C,B*N-@C, ®y, =, @C,B’@b 


@X, I B'N @y, =,B'®b 


4.1 | The Primal Simplex Algorithm of GLP 


Suppose that a basic feasible solution is also accompanied with a basis Band corresponding simplex 
table. 


I. The basic feasible solution is given by @x, =, B'®@b =, Oy, and @x,, =, @0. 


Then the grey objective function value will be @Z=, @c,-B’-@b=, @Vy. 
Il. Calculate @y,, =, @Z,-@C,, J= Ltpag hf FB Ia 1 Ay Me, 


Let @y,, = min|®y,,} 


G Isjsn 


I, IF OY,,2, @0,, then stop; the solution is optimal. 


IV. If @Y,,<, @0 and Vy. S 0 1=1,2,...,M , the problem has an unbounded solution. 


Vv. If @Y,, <6 ®0 and there is (=1,2,...,mso that Vin >0O, then determine an index r 


corresponding to a variable X,. that leaves the basis as follows: 


VI. Pivot on y,, and update the simplex tableau. Go to Step 2. 


In the following, we will look at an example to demonstrate how the above Algorithm plays out in real 
life. 


5 | Numerical Example 


In this section, by presenting example method's efficiency to solve, is evaluated. 


Example. Consider the GLP problem below. 
Maximize ®z =, ®[1,3] ®x, + [2,5] ®x, 
Subject to 

2@x, +3 Ox, <, @[5,7|, 
3 @x, + @x, <, @[3, 6], 
@x, ,@x, 80. 


Solve. The above problem is called the initial GLP problem. The standard form of the above GLP 
problem is considered as follows. 
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simplex algorithm 


Maximize ®z =, ®[1,3| @x, +@[2,5]®x, 


Subject to 
2 @x, +3 Ox, +@s, =, @[5,7], 
3 @x, + @x, +s, =, @[3,6], 
@x, ,@x,2, ©0. 
Table 2. The first simplex table of the GLP problem. 
Basic Variables @x, @x, @s, @s, Sources 
@z, -@[1,3] -@[2,5] [0,0] ef0,0] 20,0 
@s, 2 3 1 0 25,7] 
@s, 3 1 0 1 8/ 3,6] 
Table 3. The optimal simplex table of the GLP problem 
Basic Variables @x, @x, @s, @s, Sources 
e 53] eos) ofa] ope) ofa.5] 
2 1 
@x, ; 1 ; 0 2 $,2| 
7 1 
@s, 5 0 a. [2,2] 


6 | Conclusion 


In most real-world issues, we are dealing with situations where accurate data and complete information are 
not available. Therefore, LP is not accurate due to the assumption that the parameters are constant. Given 
that the theory of the GS is one of the approaches to deal with uncertainty. There is an urgent need to 
provide a suitable solution for GLP models. To date, many studies have focused on this topic, most of 
which are discussed in this article. Some of the proposed methods were not able to solve the GLP models 
in different conditions, while others did not provide the appropriate optimal solution. In general, the GLP 
problems transformed into one or a series of the classical LP problems and then obtained an optimal 
solution. As shown in the paper, an algorithm like well-known primal simplex algorithm for solving GLP 
problems is proposed, where it can solve the problem without converting to the classical LP problems. In 
future research, researchers could suggest a more accurate way to solve the GLP problem. 
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